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We study the effects of friction on the chiral symmetry restoration which may take place tem- 
porarily in high energy heavy ion collisions. The equations of motion with friction are introduced 
to describe the time evolution of the chiral condensates within the framework of the linear a model. 
Four types of friction are used to study how the result is sensitive to the choice of the friction. For 
the thermalization stage, the time dependent temperature is parameterized so as to simulate the 
result of the parton-cascade model. It is parameterized according to the one dimensional scaling 
hydrodynamics for the subsequent cooling stage. The time development of the condensates and 
the entropy production due to friction are calculated numerically. The time interval in which the 
chiral symmetry is restored approximately is investigated in detail for four types of friction. It is 
found that; (i) the maximum temperature must be high enough (not lower than 230 MeV) and the 
friction must be strong enough in order that the chiral symmetry restoration lasts for a long time 
(not shorter than 3fm/c); (ii) the ratio of time interval in which chiral symmetry is restored, to the 
time interval in which the temperature is higher than the critical temperature is typically 0.5 when 
the friction is strong enough; and (iii) the entropy due to the friction is mainly produced in the early 
stage of the cooling. The effect of freezeout is discussed briefly. 
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I. INTRODUCTION 

A new phase of matter called quark-gluon-plasma (QGP) is expected to be produced at high energy heavy ion 
collisions. The chiral symmetry may be restored in the new phase. It is usually claimed that chiral symmetry 
is restored, i.e. the chiral condensates vanish when the temperature T of the system is higher than the critical 
temperature Tc- However, this is true only when the system is always near the equilibrium. The effective potential 
may be chirally symmetric if T > Tc. However, the condition T > Tc does not necessarily guarantee the vanishing 
of chiral condensates in a time dependent process because of a finite relaxation time. In our previous work |l) where 
friction is not taken into account, it was found that the maximum temperature must be high enough for the chiral 
condensates to vanish temporarily and furthermore the chiral condensates cannot stay near the origin of the chiral 
space for a long time in general. The reason is as follows. For T > Tc, the time derivative of the chiral condensates 
becomes maximum at the bottom of the single-well potential if there is no friction. Then, it is difficult for the 
condensates to stay near the minimum of the potential for a long time. This implies that the chiral symmetry 
does not restore for a long time even when T becomes much higher than Tc- It may, however, be possible that the 
condensates stay near the minimum of the potential if there is friction of sufficient magnitude Therefore, it is 

important to study the effect of the friction on the motion of the condensates in the chiral space. Here, we consider 
that the vanishing or very small chiral condensates are the signal of the chiral symmetry restoration. 

The chiral symmetry restoration is directly related to the interesting phenomena called disoriented chiral condensates 
(DCC) H-^. In the conventional scenarios of DCC formation, it is usually assumed that the condensates roll down 
to the minimum of the temperature dependent effective potential from the top of the hill of the potential. In such 
scenarios, it is implicitly assumed that the chiral symmetry restoration takes place for a sufficiently long time. Then, 
it is obvious that the investigation of the chiral symmetry restoration in a dynamical process is important also for 
studying DCC formation. 

In this paper, we investigate whether the chiral condensates can stay near the origin of the chiral space for a 
long time when the equation of motion has a friction term using the linear a model with the massless free particle 
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approximation (MFPA) The conditions of high energy heavy ion coUisions are taken into account through an 
appropriate initial condition and a time dependent temperature. 

This paper is organized as follows. In the next section, we derive the equations of motion for the chiral condensates. 
Various types of friction are introduced phenomenologically to study its role played in the chiral symmetry restoration 



in heavy ion collision processes. An expression for the entropy production due to the friction is derived. In section II], 
the equations of motion arc solved numerically for an adequate initial condition and a time dependent temperature. 
The time interval in which the chiral symmetry is restored approximately is studied in detail. Entropy production is 
also studied. Section is devoted to conclusions and discussions. 



II. THE LINEAR a MODEL AND THE EQUATION OF MOTION WITH FRICTION 

A. Equation of motion 



The linear a model is a useful tool to investigate the time evolution of the chiral condensates. The Lagrangian is 

/: = ia^0a'>-^(02_^;2)%i7a, (i) 

where = (tr, tt) and = (0o, (pi,(j)2:4'3) = '^f=o^ 4'']- '^'^^ equation of motion for the field (pi is 

□ 0, + A _ „2) 0,^ _ HSr^o = 0. (2) 

The field (pi is divided into two parts, the condensate <I>i = {(pi) which is the expectation value of the field (pi and the 
fluctuation (pi — (pi — <I>i. Then, we obtain 



□ + A ($2 _ _ H5,^Q + + m:fp, + A 



^2 + 2^$,^, 



+A 



= 0, 



(3) 



where = 2A$f + A (<i>2 



In eq.(^, we take the normal ordering with respect to the vacuum for which the expectation value of pi is $i. 
Assuming that the fluctuation fields (pi are in the thermal equilibrium, we obtain the following equation of motion: 



A - 1,2) _ H5,o + A 



-0. 



(4) 



Note that the thermal average of the third order term in (p is vanishing under the free particle approximation. On 



the other hand, the thermal average of the 



is given by 
dk 



{2TT)^UJk,i 



exp 



1 



(5) 



where uJk,i = yk"^ + m| . It is equal to T'^/12 for rrii = 0. As a result, the equation of motion of condensates under 
MFPA S turns out to be 



One can see from eq.(@) that the potential for the condensates is given by 



T/(<i>,$o,r)-^ ($ 



We define the temperature dependent mass squared by the second derivative of V{^, 3>o, T): 



[M,($,$„r)] 



T2 

= A ( $2 - -y^ + — + 2$ 



(6) 



(7) 



(8) 



2 



B. Choice of friction 



As already mentioned in SecJ|, the aim of this paper is to investigate the friction effects on the chiral symmetry 
restoration in high energy heavy ion coUisions. The system which consists of the chiral condensates (soft modes) will 
be far from a (local) thermal equilibrium in such a time dependent process even when the environment which consists 
of the fluctuation fields (hard modes) is approximately in the (local) thermal equilibrium. In such a case, a friction 
term appears in general as a result of the dissipation- fluctuation theorem Thus we use a phenomenological 

equation of motion with a friction term: 

+ xU' + ^^ vA - + % ^ = 0- (9) 

Here, r/^'s are the friction coefficients. In general, the friction coefficients depend on both T and $i. In this paper, 
we assume that rjij = rjSij. For the functional form of 77 which may provide either T and/or $i dependence, we use 
the following form which was obtained for the chirally symmetric vacuum |2,H] (See also, p^): 



9 X'T' 

V 



167r3 



/5p(l-e-'"/^) =F(m,r), (10) 



where fsp{x) = — d.t^^^ is the Spence function and m is the temperature dependent mass. In this paper, we 
consider the following four cases as the choice of 77: 
case A : 

1] = constant. (Ha-) 

i.e. T] is taken for simplicity as a constant which is independent of both T and $i 
case B : 

77 = F(m,,r), (lib) 

where m,r is the pion mass at zero temperature. This choice may serve as a case in which rj has the simplest nontrivial 
T-dependence. This case may be realized in the rapid cooling case, 
case C : 

ry = F(Af3(/.,0,r),r). (11c) 

This is another case where rj has a simple T-dependence. In this case, we ignore the change of $ in the $- and T- 
dependent mass defined in eq.(^, i.e., we take $0 = /tt and — for i = 1, 2, 3 in order that M^{fT^, 0, T) coincides 
with the pion mass at T = 0. 
case D : 

r; = F(M3(<i>.,0,T),T), (Ud) 

where is the condensation value of the temperature dependent vacuum which is determined by the extremum 
condition 



dV 



= 0. (lie) 



In this case, the mass M3 depends on both T and $^ while the latter is a function of T. Note that only case A has 
been considered by previous authors in numerical analyses \%. 

The magnitude of the friction may be estimated from eq.(ry) for the symmetric phase. Biro and Greiner suggested 
that rj = 2.2fm~^ at T ~ Q|. Rischke calculated rj for the broken phase ||^. It was also estimated by Yabu et al. 
using the Caldeira-Legegett method. Their result is 77 = O-TttItt in the uniform case. For the case A, we use rj — 2fm 
as a typical value and consider also some other values for comparison. Considering the theoretical ambiguity in the 
magnitude of 77, we introduce the overall multiplicative constant C,, for cases B,C and D. The friction coefficient rj in 
these three cases are replaced with C^jy with various C,, in the following numerical calculations. 
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C. Time dependence of temperature 



We consider the time development of chiral condensates in the entire stage of heavy ion colHsions at high energies. 
We concentrate on the central rapidity region where the temperature will become highest. A process which takes place 
in this region will be approximately invariant under the Lorentz boost along the collision axis. Then a convenient 
variable is the proper time r. The condensates which are formed in the central region will evolve under the influence 
of the environment which consists of many number of partons (quarks, antiquarks and gluons) or hadrons. The effect 
may be represented by the temperature T which depends on the proper time r provided the environment is in the 
local thermal equilibrium . We use the temperature as a convenient parameter even for the very early stage of the 
collision process where the environment may be far from the thermal equilibrium. This is because we do not know 
any other reliable method to describe the effect of the environment in the pre-equilibrium stage. 

It is conceivable that hard collisions among the incident partons play a crucial role to thermalize the system in the 
beginning of the entire process. We then borrow the result of the parton cascade model [^3| in order to parameterize 
the T-dependence of T for the early themalization stage including the very early stage of the collision. It is possible 
to estimate the r - dependent temperature T(t) up to r = when the temperature reaches the maximum Tm using 
the energy densities calculated in with the equation of state for the free quark-gluon gas. For example, one has 
T,n ~ 620 MeV and ~ 880 MeV for Au-Au collisions at RHIC and LHC energies, respectively. The proper time Tm 
is typically some 1 fm. After this time, the system will expand mainly along the collision axis. The one-dimensional 
scaling hydrodynamics [l4| will be a good approximation for the cooling stage. Then the temperature becomes a 
function of r only, T = T{t) cx r^^/^ for r > r™. To summarize, the temperature is parameterized as follows as a 
function of the scaled proper time x — r/Tm for the entire stage of the collision [|| : 



(12) 



The temperature parameterized in eq.(p^) stays finite for a long time. However, the interactions between the 
environment and the condensates will vanish at a certain temperature called the freezeout temperature Tf. Therefore, 
it is reasonable to consider that T suddenly becomes zero when it reaches Tf. Then, the time dependence of the 
temperature may be parameterized as 



T{x)^T,{^^ e{l-x)+T^x-^/H 



[x — l)0{xf — x), 



(13) 



where Xf = {Tm/Tf) and r/ = x/t,„ being the freezeout time. To summarize, the temperature decreases from Tm 
down to Tf and then suddenly becomes zero. It is supposed that the friction vanishes after Xf because there will be 
no energy dissipation from the condensates into the thermal environment after the freezeout. In this situation, the 
initial condition of the quench scenario for the formation of DCC may be realized at a; = x^. On the other hand, 
a whole process may exhibit a characteristic feature of the annealing scenarioj7|g] if Tf is much lower than Tc and if 
the one dimensional expansion lasts for a long time. Equation of motion, eq.(0)^ is now rewritten as 



ld_ ^ 
dr^ T dr 



T\x) 



dr 



0, 



(14) 



where we assume that the system is homogeneous along the transverse direction. In the next section, eq.(p^ will be 
solved numerically with T{x) given by eq.(12) or cq.(0). 



D. Entropy production 

Consider the following equation of motion of a particle with mass m in order to find the expression for the entropy 
production: 

m^^FiX) + GiX), (15) 

where F{X) represents a conservative force and G{X) a non-conservative one. The work W which the particle does 
by G{X) is 
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W = -J^dXG{X), (16) 
where c is the path of the particle . If the work W is fully converted into heat, the produced entropy is given by 

In the present case, X corresponds to and G{X) to — riij{^)d^j/dt in eq.(||). Then, The increase in the entropy 
density may be given by 

where t has been replaced with r. The above expression is in agreement with eq.(4.6) in ref. pC| ] when there is only 
one kind of field. The entropy density production per unit proper time is given by 

dsi _ riij ( d^i \ ( 
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dT ^ T{X) \dT J \dT r ^ ' 



III. NUMERICAL RESULTS 
A. Choice of parameters 

We take the parameters of the linear a model as follows in the following numerical calculations: A = 20, v — 87.4 
and = 119 MeV. This set of the parameters generates the pion mass 135 MeV, the sigma mass 600 MeV and 

the pion decay constant f-^ = 92.5 MeV. The initial condition of the equation of motion, eq.([l^), is taken as 

($0,$) = (/.,0), (20a) 
7rr2 7rT2 T:Tr 



d$o d^i d^2 d^3 
dr ' dr ' dr ' dr 



30' V30' V30' 730 



(20b) 



where Tj is the initial temperature. 

Four types of friction given in eqs. ( 11a ) ~ ( |llcj ) are used in the following calculations. The T-dependence of 
the friction in cases B,C and D is displayed in Figs.|l^ and The T-dependence in the case D comes from the 
T-dependence of the pion mass. The pion mass is almost constant below Tc and increases linearly above T^. This 
behavior of the pion mass below Tc is similar to that in case B. The behavior above Tc is asymptotically identical to 
that in case C. Moreover, the magnitude of the friction is determined mainly by the factor T^/m(T). As a result, 
the T-dependence of the friction in case D is quite similar to that in case B at low temperatures and it is similar to 
that in case C at high temperatures. The T-dependence of the friction in case D has a peak near the Tc because of 
T-dependence of the pion mass. Since the T-dependence of the friction comes from the T-dependence of the pion 
mass, that in case B is proportional to T^, while that in case C and D is proportional to T at high temperatures. 



B. Time development of chiral condensates 

We have solved eq (^4|) numerically for various types of friction. The result for case A with the maximum temper- 
ature T,n = 20GMeV, the initial temperature Ti — IMeV , the time = Ifm and various 77 is shown in Figj^. It is 
found that a) the motion of the condensates for a; < 1 is insensitive to the value of rj , and b) the condensates cannot 
stay for a sufhciently long time in the center of the chiral space. As already found in 0, a characteristic damped 
oscillation appears in the no friction case (77 = 0). Contrary to this, the oscillation diminishes rapidly when 77 = 1 
and 2 fm^^. The condensates ($0) in these two cases pass near the center of the oscillation in the 77 = case. 

The time development of condensates in cases B,C and D with T^ = 250 MeV, Ti —1 MeV and t,„ =lfm is shown 
m FigJ|. There are characteristic oscillations in cases C and D, while there is not in case B. The oscillation in case D 
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is weaker than that in case C. These different behaviors reflect the different T-dependence of the friction in the three 
cases. 

The Tm dependence of the condensate $o in case A is shown in Fig.^, where the parameters are taken as = 1 
MeV, Tm = 1 fm and rj = 2fm^^. Higher the maximum temperature is, faster it moves to the origin. However, it is 
difficult to be stopped by friction as the condensate moves rapidly. The different behavior for a: > 10 is due to the 
difference of T. Since the behavior of T is determined by the scaling property, T depends on only if t„i is fixed. 
Higher T is, smaller the minimum of the potential (eq.(^) is. As a result, the condensate is small for high Tm- 

The sensitivity of the time development to the initial temperatures Ti is shown in Fig.^. It is found that the Ti 
dependence is weak. 

The proper time Tm is 1 fm or less according to the result of the parton cascade model. The condensates evolves 
in a different way for different even if other parameters, Tm,Ti and rj are fixed. Since determines the behavior 
of the temperature in the entire stage, we show typical dependence in Fig.|| for Tm = 250MeV, Ti — IMeV and 
77 = 2fm~^. For x < 10, the time development in x is faster for larger t„i- On the other hand, the behavior for a: > 14 
is almost independent of t„i- That is, the behavior at large x scales in x. The motion at large x is determined solely 
by T for given Tm, T and rj. 

The effect of freezeout on the time development of condensates in case D is demonstrated in Fig.|^. The curves (a) 
and (/3) represent the results obtained with the time dependent temperature ( jl^) and ([l3|), respectively. Of course, 
the effect of freezeout is seen only for x > Xf. For comparison, T{x) given by eq.(|l3|) is also shown in Fig.M. The center 
of the oscillation of sigma condensate in the curve (/3) is almost The amplitude of the oscillation is determined 
by the magnitude of the condensate at the freezeout time. 



C. Chiral symmetry restoration time 

One may claim that the chiral symmetry is restored approximately if every | $i | becomes much smaller than v. 
For definiteness, we use the following condition to define the approximate restoration of chiral symmetry: 

V$2/w<0.2. (21) 



From now on, we say that chiral symmetry is restored if the condition ( |21| ) is satisfied. 

We take Ti = 1 MeV and r„j = 1 fm in the following calculations. Time interval in which chiral symmetry is 
approximately restored is calculated for the four cases with various Tm- The result for case A is shown in Fig.|[ It 
is found that the condition Tm > 230 MeV, rj > Ifm^^ must be fulfilled in order that > 3 fm. At fixed 77, higher 
Tm is, longer becomes. This is because the condensate follows the minimum of the potential realizing the chirally 
symmetric vacuum and high temperature era lasts for a long time when Tm is high. On the other hand, for fixed 
Tm, Tr increases as a function of rj for small 77, reaches a maximum at some rj and then decreases. This behavior 
is understood as follows. If 77 is large, the condensate moves slowly and almost stops near the origin resulting in a 
long Tr- However, if 77 is too large, T decreases considerably before the condensate reaches the origin and hence, r,. 
becomes short. The same tendency is found in other cases. 



The result for the case B is shown in Fig.^ . Now Tr is plotted as a function of Ct^ introduced in subsection pB . 

The condition for having Tr > 3fm~^ is Tm > 250 MeV and ~ 0.1. In this case, large can be realized even when 
Cjj is very small because 77 becomes very large. (See fig.|^) The result for case C is shown in Fig.|l0|. The appearance 
of large for Tm = 140 MeV is due to an accidental fluctuation. Maximum temperature must be higher than 260 
MeV and must be larger than 2 in order that > 3 fm. Fig.|ll| shows the time interval in case D. The condition 

for having > 3 fm is T„ > 220 MeV and - 1.5 - 2.0. To summarize the results for the four cases, we found the 
following general tendency common to every case: (a) the Tm dependence of is approximately universal as long as 
the friction is sufficiently but not too strong; (b) Tm must be larger than 220 ~ 260 MeV in order to have > 3 fm; 
(c) approximate chiral symmetry restoration does not take place if the friction is very small. 

We define the time interval Tc in which the temperature of the system is higher than the critical temperature Tc in 
order to compare it with r^. It is given by 



Tm 

T, 



(22) 



This ratio in case A is shown in Fig.|T^. It is found that the typical value of this ratio is about 0.5. Each of the lines 
corresponds to various Tm ranging from 230 McV to 390 MeV. The ratio is approximately independent of 77 for large 
77. We have confirmed that these features are common to all the four cases. 
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D. Entropy production 



Finally, the entropy production due to the friction is calculated. The result for case D is already shown in figj^- 
For comparison, we have calculated it for case A with i] = 2fm~^, Tm = 250 MeV, — 1 MeV and Tm — 1 fni. The 
result is shown in Fig.|l^. It is found that entropy is produced mainly in the early cooling stage. Entropy production 
lasts for a relatively short time (a few fm) and hence it is not affected by the freezeout. The entropy production is 
mainly due to the motion of <i>o condensate 



IV. CONCLUSIONS AND DISCUSSIONS 



The time development of chiral condensates in high energy heavy ion collisions are calculated with the adequate 
initial conditions and frictions. The linear a model with four types of friction are used. The main results are: a) 
the characteristic damped oscillation of the condensates in no friction case is smeared by the friction; b) sufhciently 
high maximum temperature (T,„ > 220MeV) and the friction of appropriate magnitude are needed for the vacuum 
to be chiral symmetric for a sufficiently long time (r,. > 3fm); c) entropy is produced mainly in the early stage of the 
cooling. 

For each friction type, we found the following results. In case A, the constant friction which has been used by some 
authors generates the largest among four cases. However, it is not clear whether the T independence of ry is a 
good approximation or not. In case B, is long enough even for small C'ri- However, the magnitude of the friction is 
unrealistically large at high temperatures. In case C, is small for the realistic range of the friction (0.5 < Cjj < 3) 
. In case D, large is obtained for the possible range of the friction (0.5 < < 3). The peak of 7] near Tc (see 
Fig.pH]) is responsible for large in this case. 

As stated in the introduction, the conventional scenarios of the formation of QGP and DCC in which the chiral 
symmetry is supposed to be restored for a sufficiently long time are not always guaranteed. Our results suggests that 
it is important to estimate precisely the maximum temperature of the system and the friction for the soft mode. 

In general scenarios of QGP evolution, it is well-known that the entropy of the system is mostly produced in the 
thermalization stage (x < 1). However, we found that the entropy due to the friction felt by the condensates is 
produced in the early stage of the coohng. 

As the initial condition of conventional DCC formation scenarios, it is often assumed that the roll-down of chiral 
condensates starts from the top of the potential (chirally symmetric vacuum) and the initial temperature is close to 
Tc- However, this initial condition may not be always satisfied as the chiral symmetry restoration does not necessarily 
occur even when the temperature of the system exceeds Te- 
la our calculation, the condensation dependence of the friction is included in case D. However, it reflects essentially 
only the T-dependence of the vacuum. The effect of the time dependent condensates is not taken into account yet. 
Such a dependence |3| should be included to refine the analysis. On the other hand, random forces are not considered 
also in this paper because we are interested in the average behavior of the chiral condensate. They generate fluctuations 
in the motion of the condensate and hence in r^. The fluctuation-dissipation theorem implies that the random force is 
strong when the friction is strong. However, the strong random force does not always affect considerably the motion 
of the condensate because its effect is reduced by the strong friction. Even when the random force works, Tr will 
become shorter on the average because condensate will occasionally be kicked out from the chirally symmetric region 
by the strong random force. This means that higher Tm is required to realize the same on the average. There is 
another aspect in this situation. The condensate in the present case moves mainly along the a axis. It is almost at 
rest along the tt axis. Though the random force will decelerate or accelerate the condensate along the sigma direction, 
it will only accelerate the condensate along the if direction. As a result, the chiral symmetry will be more easily 
broken when the random force is included. 

Finally, we would like to make a comment on the massless free particle approximation. It is not reliable at low 
temperature. The effects of mass and the finite volume of the colliding system have to be taken into account to 
improve the theory. However, we believe that our results are correct at least qualitatively. The problems listed above 
and their effects on DCC formation will be discussed in the coming paper. 
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FIG. la. The temperature dependence of the magnitude of friction in cases B,C and D. 

FIG. lb. The temperature dependence of the magnitude of friction in cases B,C and D. Note that the scale of T is different 
from that in Fig.[la|. 

FIG. 2. Time development of chiral condensates for various magnitudes of friction with Tm = 200 MeV, Ti — IMeV and 
= Ifm in case A 



FIG. 3. Time development of chiral condensates for Tm ~ 250 MeV, Ti — IMeV and Tm = Ifm in cases B,C and D. 
FIG. 4. Time development of chiral condensates for various Tm with T = IMeV, Tm = Ifm and t] — 2fm~^ in case A. 
FIG. 5. Time development of chiral condensates for various T with Tm = 250 MeV, Tm = Ifm and 77 = 2fm^^ in case A. 
FIG. 6. Time development of chiral condensates for various Tm with Tm ~ 250 MeV, T — IMeV and rj — 2fm^^ in case A. 

FIG. 7. Time development of sigma condensate and entropy production per unit scaled time (r/r™) without (q) and with 
(/?) the freezeout in case D. The parameters are Tm = 250 MeV , Ti = 1 MeV, r™ = Ifm and Tj = 90MeV. 

FIG. 8. The time interval (r^) in which chiral symmetry restoration lasts longer than 3fm for various Tm with T = 1 MeV 
and Tm ~ i fm in case A. 
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FIG. 9. The time interval (r,.) in which chiral symmetry restoration lasts longer than 3fm for various Tm with T, = 1 MeV 
and Tm, = 1 fm in case B. it is displayed for (a) < C,, < 1 and (b) < < 3. 



FIG. 10. The time interval (xr) in which chiral symmetry restoration lasts longer than 3fm for with Ti = l MeV 

and Tm = 1 fm in case C. it is displayed for (a) 0.5 < Cr, <3 and (b) < Cn < 20. 



FIG. 11. The time interval (xr) in which chiral symmetry restoration lasts longer than 3fm for various Tm with T, = 1 MeV 
and Tm = 1 fm in case D. it is displayed for (a) < C,, < 3 and (b) < < 20. 



FIG. 12. The ratio Tt-/tc for various Tm in case A. 



FIG. 13. Entropy production per unit scaled time {r /Tm) without the freezeout with Tm = 250 MeV , Tj = 1 MeV, Tm = 
Ifm and rj = 2fm~'^ in case A. 
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